LOW TEMPERATURE VISCOSITY IN ELONGATED FERROFLUIDS 



T. Alarcon, A. Perez-Madrid and J. M. Rubf 

Departament de Fisica Fonamental 
Facultat de Fisica 
Universitat de Barcelona 
Diagonal 647, 08028 Barcelona, Spain 

We have studied the relaxation and transport properties of a ferrofluid in an elongational flow. 
These properties are influenced by the bistable nature of the potential energy. Bistability comes 
from the irrotational character of the flow together with the symmetry of the dipoles. Additionally, 
the presence of a constant magnetic field destroys the symmetry of the potential energy magnetizing 
the system. We have shown that at a moderate temperature, compared to the height of the energy 
barrier, the viscosity decreases with respect to the value it would have if the potential were stable. 
This phenomenon is known as the 'negative viscosity' effect. Thermal motion induces jumps of the 
| magnetic moment between the two stable states of the system leading to the aforementioned lowered 

dissipation effect. 
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I. INTRODUCTION 



> 

■ The existence of a 'negative' viscosity in dispersions of ferromagnetic particles in a nonpolar .solvent is a curious 
' phenomenon which has recently been discovered. It was predicted by Shliomis and Morozovu and corroborated 
experimentally by Bacri et aln. Its interest is based upon the fact that contrary to what one would expect, the 
viscosity of the dispersion as a whole diminishes when an applied magnetic field oscillates at an optimum frequency. 
As has been known for a long time, the presence of a constant magnetic field prevents free rotation of the dipoles 
which originates a resistance on the flow of the fluid which increases dissipation. On the other hand, Bacri et aln have 
shown that the increment in the viscosity due to the rotational degrees of freedom is proportional to the difference 
between the vorticity and the angular velocity of the particles in suspension. Thus, a magnetic field oscillating at a 
high enough frequency can impart an angular velocity to the particles greater than the vorticity, leading to a 'negative' 
I . viscosity. One see that the particles gain kinetic energy at the expense of the oscillating field. 

This striking phenomena is the motivation for undertaking this investigation. We have found that the effect 
described is more general than was thought, considering that it can be applied to any ferrofluid that has a frequency 
O ■ different from the inverse of the Brownian relaxation time. 

The key point in our results comes from the fact that we work with axisymmetric particles immersed in a suctioning 
current. Thus, at a low enough temperature, due to thermal agitation, the particles jump back and forth along the 
k><| \ axis of the flow, which thereby introduces an internal frequency in the system. This internal frequency plays exactly 
the same role as the frequency of the alternating magnetic field in the systems studied by Shliomis and MorozovEl and 
' Bacri et ala. 

There are examples of systems other than ferrofluids that potentially can show this phenomenology, these are 
suspensions of gravity dipolestm or suspensions of magnetotactic bacteriaB. These bacteria, some of them rod shaped, 
undergo the phenomenon of magnetotaxis. The bacteria contain iron particles which impart a magnetic moment to 
themselves. 

The paper is organized as follows; in section II we describe the system and perform an analysis of the stability of 
the equilibrium orientations of the dipoles. Section III is devoted to the derivation of the relaxation equation for the 
magnetization. In section IV we introduce fluctuations in our analysis and, by applying fluctuating hydrodynamics 
in the space of orientations, we derive the Langevin equation for the magnetic moment. In section V we compute the 
stress tensor and the viscosity tensor. Finally, in section VI, we discuss our main results. 



II. THE MAGNETIC ROTOR 



The system we want to study consists of a dilute colloidal suspension of ferromagnetic rodlike particles immersed in 
a nonpolar fluid phase. This suspension flows in an elongational flow and under the influence of a constant magnetic 
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field H oriented in the direction of the symmetry axis of the flow, which will be taken parallel to the z-axis. All of 
the particles are supposed to have the same magnetic moment 

m = m s R, (1) 

where R is the unit vector along the direction of the axis of the particle and m s is the magnetic moment strength. 
The velocity field of the flow is given by 

<7=|-r, (2) 

where (3 is the velocity gradient of the flow 

* (- 1 °\ 

/? = /? -1 (3) 
\ 2 / 

with f3 being the rate of elongation. 

A suspended particle in the carrier fluid experiences the hydrodynamic torque given by 

f h = -t-(£3-Si), (4) 

where the friction tensor £ is given in terms of its components £o and £i 

t=£iRR + £o(f -RR). (5) 

Here uj, ft and 1 are the angular velocity of the particle, the drag angular velocity due to the motion of the fluid and 
the unit tensor, respectively. Due to the fact that the only possible relative motion between the ends of the rod is a 
rigid rotation one has 

v+ - v- = ft x RL, (6) 

where v± are the velocities of the end points of the rod and L is its length. This expression can alternatively be 
written as 

v+-v_=f3-RL. (7) 
Thus, through a comparison of eqs. (||) and ([?]) we achieve 

ft = 7?(l/2/3 : RR), (8) 

with TZ = R X -^g being the rotational operator, 
i-i 

In addition to eq. (|4|) the dipole is acted upon by a magnetic torque 

f m = m x H. (9) 

Therefore, by taking into account the angular momentum conservation, in view of eqs. (^), (^), (j^), and (^J) we can 
write in the high friction limit 

ts=h-nu), (io) 

?0 



where the potential U, defined through 



U = -m ■ H - iftl : RR; (11) 



accounts for the different mechanisms for the rotation of the particle: the magnetic field and the external flow. In 
polar spherical coordinates this energy can be rewritten as 
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U = -m s Hcosd+ ^-(1 - 3 cos 2 (9), (12) 

8 being the angle between the axis of the particle and the z-axis. 

The equilibrium states of the system can be identified through the condition = 0, i.e. 

m s H sixiO + 3f /3 sin 9 cos = (13) 

whose solutions, 9 + and 8q, are 

9- = 0, (14) 

0+ = 7T, (15) 

0o=arccos (~il)- (16) 

The stability of these orientations follows from the second derivative of the potential computed in each of the solutions 

^\e_ =m s H + 3£ /3 = m s H(l + ^), (17) 

-^-\ e+ =-m s H + 3Zo(3 = m s H(^-l), (18) 

^ k_ ^r - ^- (19) 

Since m s , H, £o and f3 are always positive quantities, after examining eqs. (0)-(^9|) we conclude that provided 

H < H c , (20) 

H c = 3£ f3/m s being a critical field, 6L and 8 + are both stable, whereas Qq is unstable. We then conclude that under 
these conditions the potential of the magnetic rotor is bistable. For magnetic fields larger than H c the system becomes 
stable. 

Equivalently, we could have considered the presence of a critical value for the elongational rate, (3 C — ^ , such 
that for a fixed value of the magnetic field the system is bistable whenever the actual value of the elongational rate 
j3 overcomes (3 C . In the rest of the paper we will assume that H/H c < 1, i. e. we will remain in the bistable region, 
or in the range of large elongational rate. 

To conclude this section we will give some estimates of the critical field corresponding to situations of experimental 
accesibility. For particles of magnetite having a volume V p — 5 x 10~ 19 cm 3 and an aspect ratio e = 0.1 one has 
m s = 2.4 x 10~ 16 G x cm 3 , and consequently H c — 58.30e. On the other hand, for particles of cobalt with volume 
V p — 2.7 x 10~ 19 cm 3 and e = 0.1 we obtain m s — 3.8 x 10 _16 G x cm 3 , and H c = \7 .lOe. In both cases we have 
assumed (3 — 10 3 s _1 . At smaller elongational rate, (3 = 10s _1 , one has H c = 5.80e for magnetite and H c = l.70e 
for cobalt. 



III. RELAXATION EQUATION 



The relaxation of the magnetic rods can be interpreted as a diffusion process through a potential barrier in orien- 
tation space. In this scenario we can apply the formalism of non-equilibrium thermodynamicsa. Let p(R, i) be the 
angular distribution function that may be viewed as a density in the space of orientations. Thus, we can define the 
chemical potential as 

fi = K B T\np+U. (21) 
The diffusion equation for p is written in the form 
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dp D d , , n , p dU dp,, , . 

Tt = ^rede^ 6 ^m + i^ (22) 

with D = ksT/^Q and where the axial symmetry of the potential ( |l2| ) has been taken into account. Equation (|2 
can be rewritten as a continuity equation 

which defines the diffusion current 

Jo = -De(- u ' K » T ^e»l KBT . (24) 
09 

If the height of the potential barrier is large enough as compared to thermal energy KbT, we can suppose that 
equilibrium is reached independently on each side of the barrier. Thus, the chemical potential can be inferred as being 

p(R, t) = m(0-)6(0 o - 0) + p{9+)®{9 - d ), (25) 

where 0(9) is the unit step function. The system is allowed to obtain the global equilibrium due to the presence of a 
quasi-stationary current J(i) , assumed uniform. 

J e sm9 = J(t){G(0 - 9-) + Q{9 - 9+)}. (26) 

This adiabatic hypothesis is one of the essential points in the present development. By using (|2l]) and (p5|), one has 

P (R,t) = p-(t)el-( u - u -y K * T »e(6 -9)+ p + (t)e^ u - u +V KBT Ve(6 - 9 ), (27) 

where p± = p(9±,t), and U± = U(9±), are the densities and potential energies in the two stable states. In order to 
obtain an expression for J(t), we substitute ( p4] ) into (|2^) to arrive at 

P U/K B T fj 

J{t) -^9~ {Q{B 9 - ] + 0(6 9+)} = ~ D d9 e ■ (28) 

Integrating now over 9 and taking into account that, due to the height of the barriet. the main contribution to these 
integrals is around the maximum of the potential 9q, one has the law of mass actiontO 

J{t) = K B l{l-e A ' KBT ) 1 (29) 

which is a nonlinear phenomenological relationship between the quasi-stationary current J(t) and the affinity A = 
/(+ — where p± = p(9±). The phenomenological coefficient I is given by 

I = ^^P-i^^) 112 ^-^ 1 ^, (30) 
Kb zitKb-L 

where [7 ' = §^\e and U Q = U(9 ). 

We are now prepared to proceed to the deduction of the relaxation equations. To this end, we define the following 
populations: 

N + = / p(R,t)dR= / 2tt sin 9 p(R,t)d9, (31) 



and 



0=o Je 



2-n r e a 



N- = I I p(R,t)dR= [ ° 2tt sin 6p(R, t)d9, (32) 
J0=o Je=o Jo 

related by the normalization condition 

N = N + +N_. (33) 
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The kinetic equations for N + and iV_ follow after differentiating the eqs. ( |3l| ) and ( p2\j and employing eq. (23). Thus, 
we have 



dN, 



dt dt 

where, consistently with the adiabatic approximation 

2-kKbT 



= 2nJ(t), 



U_ 



(34) 



(35) 



N+ = —^~P+, (36) 
in which U ± = ^r\e ± - Using equations ( p9| ) and (pi[)-(p6[) we can derive the rate equations for the two populations 

where the rate constants are given by 

= £sm% JU^ f {Uo _ u )/KbT 

+ 2TvK B T y 2irK B T J ' V ' 



and 



= Dsm^o ( J^L )l/2 u" e -(u -u + )/K B T_ (39) 

+ 2ttK b T { 2ttK b T j + y ' 

We can now define the magnetic moment parallel to the direction of the applied magnetic field as 

N + - N - 

m = m s — . (40) 



By differentiating equation (40) and employing equation ( p7| ) we obtain the relaxation equation for m 

dm 11 

— - = m H m s , (41) 

dt t a 



where 

1 



K+- + K_ + 

is the relaxation time, and 

1 



K+- - K_ + 



(42) 



(43) 



accounts for the asymmetry of the potential. An equation similar to ( [i"l| ) was postulated by Shliomis in the context 
of ferrohydrodynamicsQ. Our approach to deduce eq. (41) has been based upon mesoscopic arguments and as we will 
see in the next section it provides the natural way of introducing fluctuations in the scheme. 

It is useful for our purposes to introduce the nondimensional variables x = and /i = ™ 3 ^ , in terms of which the 
relaxation time is written as 
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(27r) 3 /2 r A*,l 

- 1) exp{ M (l - |)} + n(x + 1) exp{- M (l + I)}]" 1 , (44) 

where x must be greater than the unity. The inverse of r gives us a characteristic frequency proper to the system, 
the jump frequency between the two stable states of the potential energy. This frequency is given by 

w = jD (^)3/2 (1 _ _L) e -¥ e -/*/*[e^ + e"" + I( e -" - e")l. (45) 
Z7r or a; 

In Figure 1 we can observe the behavior of u> when one varies the elongational rate. The way in which it depends 
on the elongational rate and the existence of two time scales in our problem, Brownian D~ x and r, are the key points 
in understanding the dynamical mechanism which leads to the results we will obtain in section 5. 



IV. DYNAMIC OF THE FLUCTUATIONS OF THE MAGNETIC MOMENT 



Thermal motion of the dipoles inside the carrier fluid produces fluctuations in the population of the minima of the 
potential. This fact manifests itself on a mesoscopic level through fluctuations of the magnetic moment. These can 
be taken into account by adding a random currents to equation (p9f), 

J(t)=K B l(l-e A / K * T ) + — £(t), (46) 

m s 

where D is given by 

D = 2m 2 s k B l, (47) 

I is the equilibrium value of the phenomenological coefficient I, which according to ( p0|) should be proportional to p_ 
computed at equilibrium, pi 9 . By applying the detailed balance principle to eq. (p7|) and using the normalization 
relation eq. (^) one can easily achieve 

eq _ ( N \ Ulu'j 

p - \2<Kk B T) c/; + c/ : CX p { i^± } - l48j 

Additionally, in eq. ( f46| ) £(t) is a Gaussian white noise stochastic process of zero mean and correlation function. 

mm) = s(t - f). (49) 

From equations ([34]), (|4(]) and ( (46j ) we obtain the Langevin equation for m 

^ = --m + -m s + y/D^t)- (50) 
at t a 

Notice that coefficient D is the input noise strength corresponding to the stochastic process m. This coefficient can 
by explicitly written as 



and in terms of the nondimensional parameters x and fi 

fl = (4^^) 1/2 ( 1 - - ^ (l+^U+V-l)^ CXP{ 4 ( ^ + 3 )}6 ^ {/i(1 " 3 )} ' (52) 

Finally, in reference to eq. (pOj), by performing the change 
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7" 

rh = m m s , (53) 

a 

this equation becomes 

^=--m+VDm, (54) 

at t 

i.e., rh is an Ornstein-Uhlenbeck processi. As is well known, the stationary distribution for such a process is 

I TV? 

p(m) = ^==exp{ =-}, (55) 

VittD tD 

which will be used in the next section to compute the viscosity. 

V. THE VISCOSITY TENSOR 

In order to calculate the viscosity tensor, we first have to compute the stress tensor. This quantity has two 
contributions, one that comes from the solvent, and the other due to the presence of the particles. The latter is 
writtenEj as 

a = nK B T{3(RR) - f ) + <i/3 : (RRRR) - nK B Tii{RH ■ (f-RR)), (56) 

where n is the concentration of suspended particles and L is their longitude. Thus, eq. (^6|) give us first order 
contributions to the viscosity tensor. The moments in eq. ([56]) will be calculated by using the stationary distribution 
([55|). It is important to keep in mind that although we are using a stationary distribution, the dynamical effects are 
taken into account through their dependence on the relaxation time r. 

We will illustrate the behavior of the viscosity by explicitly computing the parallel viscosity defined through 

m\\ = a -f (57) 

with = H ■ a ■ H , that from eq. ( |56|) is written as 

am, = nK B T(3(R {l R {l ) - 1) + ^((R^) - (R^Rfy)) - nK B Tn({^) - (A^&q)). (58) 

Here we have taken into account that /3 : RR = (3(1 — 3i?|), and i?m| = R ■ HH. In the appendix we summarize the 
result of the computation of all the moments that appear in the expression for <7||m. Making use of these results, we 
find that 



18-— 3 - 3(-) 4 ] - nK B Tfx[- - - 3 - 3~— 5 . (59) 
a Znig a a a a 2m s 

Thus, wc finally obtain 
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1 t tD t tD 



T 



3n£o #M a 2to^ a 2777,2 q. 



1 £i r r-D 1 o r t_D n/ rD 2 . r . 2 „ , t , 4 



18-— - 9(— )' + (-)^ - 3(-) 4 ]. (60) 



3 £o 2rTig a 2m 2 , 2m 2 

3n£n 



In figures 2 and 3 we have plotted, the nondimensional quantity fjp ar = g— t in terms of x and w, respectively, for 



particles with an aspect ratio e = 0.1, for which one has i|i = 0.18050. 

One can observe that this viscosity becomes negative for large shear .rates and saturates at a positive value. Similar 
results have been obtained by Shliomis and Morozovtl and Bacri et alB with ferrofluids under an oscillating magnetic 
field. Negative viscosities in these cases were due to the existence of two characteristic time scales in the system that 
enter in competition. One of them is related to the frequency of the alternating magnetic field and the other is related 
to the vorticity of the fluid. In our case, the dipoles relax to the equilibrium orientations in a time scale D , which 
is shorter than the relaxation time r associated with the diffusion through the potential barrier. The rods are then 
constantly jumping with a frequency u> = r _1 acquiring a net angular velocity different from zero, whose expression is 
given in the appendix. The existence of this angular velocity then shows the conversion of thermal energy into kinetic 
energy for rotation. This explains why, in our case, the viscosity diminishes. 

In looking at fig. 3 one first observes a histeresys cycle described by rj par when u varies. During the first stage, when 
the actual value of u> is between zero and its maximum value, rj par , decreases from a positive value to a minimum 
negative value, as occurs in the case analyzed by Shliomis and Morozova. These results are in good agreement 
with those found by them. The difference between our approach and the one by Shliomis and Morozovu is that the 
externally imposed frequency can be as large as you want. In our case the frequency of relaxation, the result of an 
internal dynamical mechanism, does not grow freely. In the second stage; a; is a decreasing function of the shear rate. 
This fact is responsible for the particular behavior of rj par observed. Note that, as can be seen in Figure 2, r\ par does 
not saturate to a negative value when (3 — > oo. It tends very slowly to a positive value_as slowly as u> tends to zero. 
Our results also agree with the analysis of the phenomenon carried out by Rosensweigij as well. 

VI. CONCLUSIONS 

In this paper we have shown that in a ferrofluid made from rodlike particles under a constant magnetic field and 
in a elongational flow, rotations of the particles lead to aon-monotonous behavior of the viscosity. This fact allows us 
to generalize the phenomenon discovered by Bacri etaln. 

As a consequence of the orientating effect of the flow, the particles rotate with a drag angular velocity 51 that 
eliminates the degeneracy of the direction of their actual angular velocity in a volume element of the ferrofluid (see 
eq. ([To|)). Thus, all the contributions to dissipation in a volume element add up constructively. 

Moreover, because the flow is elongational, we can write a potential energy of orientation related to it. This 
potential is bistable and the addition of a magnetic field breaks up its symmetry, thus magnetizing the system. On 
the other hand, thermal motion causes jumps between the two stable states of the potential with a certain frequency. 
Consequently, the particles acquire a kinetic energy for rotation at the expense of thermal energy. This fact eliminates 
the rigidity that the system would have if no motion were present, i.e. if the potential time independent were stable. 
In fact, this is the way that we interpret the 'negative viscosity' effect predicted for our system. 

We have derived the kinetic equations for the population of the minima of the potential and from it the relaxation 
equation for the magnetic moment. Likewise, by assuming fluctuations of the density of states in the orientation space, 
we have formulated the Langevin equation for the magnetic moment. This equation describes a Ornstein-Uhlembeck 
process whose moments are well known. The computation of the first four moments allows us to derive the viscosity. 

It should be emphasized that the viscosity is computed in the limit of high energy barxiex, i.e. for strong rate of 
elongation. This is the opposite situation to the case studied in the two previous papersEd'Ei where we covered the 
weak flow regime, showing an increase of the viscosity. Analogously, modifying the magnetic field without altering 
the rate of elongation we achieve the same effect, that is, to vary the viscosity. Thus, a possible application of our 
results is in adaptive dampers. 

Through this behavior of the viscosity, emerge the macroscopic consecuences of the dynamical bifurcation with 
exchange of stability that our system experiences. 
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Likewise, the effect that we are studying is only possible in dilute solutions, where each rod can rotate freely without 
interference by others. Of course in neglecting hydrodynamic interactions among the particles we assume a extremely 
dilute ferrofluid. Anyway, this is the first step in our study of nonlinear or hysteresis effects in the rheology of 
fcrrofluids suspensions. 

The following step will be to assume higher concentrations. Nonetheless, to include hydrodynamic and exclude 
volume interactions in the dynamics of an assembly of rod particles is difficult. So that, we will model the elongated 
dipoles by means of rigid dumbbells, such models are well studied in the field of polymeric liquids. 

In order to increase the modulation effect of the magnetic field or of the elongation rate on the viscosity, another 
possibility we are thinking of consists in adding nonmagnetic spheres at high concentration to the ferrofluid. 
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APPENDIX A: 

To begin with, in this appendix we give the moments which must be calculated in order to compute the stress 
tensor. 

= - (Al) 

a 

(R^ =^=r+ - = wr+ - A2 
a Zm^ a 

{RRR)=3 lRL + { L^ (A3) 
2mj a a 

{R\\R\\R\\R\\) = 3(^) 2 + 6(-) 2 ^ + (-) 4 . (A4) 
Zmi a Zmi a 



Finally, according to eqs. (|8j), (|10|) we can write 



lj - d = —rh x H. (A5) 



Now we define the root mean square angular velocity oj rm = y ((to — ft) 2 }, that in view of eq. ( |A5|) is given by 

f tD /t\ 2 1 1/2 
u r m = nD\l-— s -(-) . (A6) 
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FIGURE CAPTIONS 

• Figure 1.- Nondimensional frequency versus the scaled elongational rate x, for /i — 1. 

• Figure 2.- Normalized viscosity versus the scaled elongational rate x, for fj, — 1. 

• Figure 3.- Normalized viscosity versus the nondimensional jump frequency for ji = 1. 
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